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ABSTRACT 

We present integral field spectroscopic observations of the quadruple-lensed QSO SDSS 1004+4112 
taken with the fiber system INTEGRAL at the William Herschel Telescope on 2004 January 19. In 
May 2003 a blueward enhancement in the high ionization lines of SDSS 1004+41 12A was detected 
and then faded. Our observations are the first to note a second event of similar characteristics less 
than one year after. Although initially attributed to microlensing, the resemblance among the spectra 
of both events and the absence of microlensing-induced changes in the continuum of component A 
are puzzling. The lack of a convincing explanation under the microlensing or intrinsic variability 
hypotheses makes the observed enhancements particularly relevant, calling for close monitoring of 
this object. 

Subject headings: gravitational lensing, quasars: emission lines, quasars: individual (SDSS 1004+4112) 



1. INTRODUCTION 

The large separation (14.62 arcsec) quadr uple gravita- 
tionally lensed qua sar w as first identified bv llnada et alJ 
( 2003) (hereinafter in the Sloan Digital Sky Survey. 
The derived redshift is z = 1.732 and the lensing complex 
consists of a cluster of galaxies at z = 0.68. 

The components of a multiple lensed system can 
show uncorrelated variability due to the gravita- 
tional (de)magnification produced by compact deflectors 
aligned with the source. This p henomenon, microlensing 
ijChang fc Refsda]|l"l979l I1984|) . is routinely observed in 
the continuum emission of multiple imaged QSOs, and 
it is considered a powerful tool for ma pping the QSO 
continuum source fYonehara et al.l fl999) . 

According to our current understanding of the struct- 
ure of active galactic nuclei (AGNs), the continuum 
source of a QSO is surrounded by a wider region, where 
the broad emission lines (BELs) originate. Early ver- 
sions of the standard model of AGNs assumed a too 
large size for the broad line region (B LR) (~1 pc) to 
be microlensed by solar-mass objects (Ncmiroff Il988t 
ISchneider fc Wambsganssl lT990). However, these early 
estimates were reduced by two orders of magnitude by 
the new measurement s of the BLR size based in the 
rever beration method ijWandel et al J 119991 IKaspi et al.1 
2000). Considering the new size estimates, a variety 
of theoretical models for the kinematic structure of the 
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BLR llPopovic. Mediavilla fe Munoz H200lUAbaias et alJ 
2002] ILewis fc Thatal2 004^ have shown that microlensing 
by stellar sized objects can produce significant amplifi- 
cations of the BEL of multiple imaged QSOs, especially 
i n the high ionization lin es. 

iRichards et alJ l)2004bl) (hereinafter IR04T ) found excess 
in the blue wings of several high ionization lines of SDSS 
1004+4112A, the highest magnified lensed image, rela- 
tive to the same lines in the other images. The excess 
persisted for at least 28 days and then faded. No such 
enhancements in the blue wings were seen in the low ion- 
ization lines. These results are in agreement with current 
ideas about the str atification of the BLR acc ording to 
the ionization level l|Peterson fc Wande]ll2000T) and sug- 
gest the possibility of scanning the BLR with high spatial 
resolution. There is no evidence of continuum microlens- 
ing in the SDSS 1004+4112 event (however see lOta et ail 
2006). This lack of correlation can b e explained for a 
particular event ({Lewis fc Ibatall2004f) . but the appar- 
ent absence of continuum amplifica tion in all the re- 
ported cases of BEL micr olensing l)R04llChartas et aD 
2002;, l200l IDai et aU feOOai is intriguing. 

In this article we present new spectroscopic measure- 
ments of SDSS 1004+4112 that show a recurrence of the 
blueward magnification of the high ionization line wings 
of image A. 

2. OBSERVATIONS AND DATA REDUCTION 

SDSS 1004+4112 was observed on 2004 January 1 9 
with the fiber system INTEGRAL ijArribas et al.lll998|) . 
at the 4.2 m William Herschel Telescope (WHT, Roque 
de los Muchachos Observatory, La Palma, Spain). Two 
sets of three 1800 s exposures each were taken co- 
rresponding to two pointings separated by few arc se- 
conds. We use the standard fiber bundle #3 (SB3) which 
consists of 115 object + 20 sky fibers, each one sub- 
tending 2.70" in diameter and covering a sky area of 
33.6" x 29.4". A 300mm" 1 diffraction grating was used, 
giving a spectral coverage of ^6000 A centered at 5700 
A and a spectral resolution of 19 A, with a spectral dis- 
persion of 6.2 A pix -1 . The observations were performed 
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TABLE 1 

Summary of SDSS 1004+4112 observations 



Date 



Telescope Ref. 



Components 



2003 May 31 
2003 Nov 21 
2003 Nov 30 
2003 Die 1 

2003 Dec 22 

2004 Jan 19 
2004 Mar 26 





174 
183 
184 
205 
233 
300 



Keck 

ARC C 

ARC 

ARC 

ARC 

WHT 

ARC 



R04 
R04 
R04 
R04 
R04 

R2 b 



ABCD 

AB 

AB 

AB 

AB 

ABCD 
AB 



a Elapscd days from first observation. 
b See lR.ichards et alJ I20n4afl 

c Astrophysical Research Consortium (ARC) Telescope 



under medium seeing (1.5"-2") conditions. The data re- 
duction procedures included cosmic ray rejection, aper- 
ture tracing, spectra extraction, wavelength calibration, 
throughput correction, sky subtraction, atmospheric ex- 
tinction correction and flux calibration. Since the vast 
majority of the flux is concentrated in one or two fibers, 
the spectrum has been obtained directly by combining 
their fluxes. 

The absolute flux calibration must be treated with 
caution since the calibration star position within the fiber 
face could affect the amount of flux lost in the intcrfi- 
ber regions. The goodness of our relative photometry 
is shown by the comparison with the photometry of 103 
(Figures and |1J. 

Apart from WHT observation we have use data from 
the Keck and ARC telescopes (See Tabled, which has 
been kindly provided by Gordon T. Richards. 




1400 1600 
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Fig. 1. — Comparison between spectra of components A (black 
in both panels), B (upper panel), and C (lower panel) for SDSS 
1004+4112 at the WHT. Spectra B and C have been scaled to nor- 
malize its continuum to that of A in the 1700-1800 A wavelength 
range. 



3. RESULTS 

3.1. Comparison of A, B, and C spectra 

Inspection of Figure d reveals that: i) component A 
shows enhancements relative to components B and C in 
the blue wings of the high ionization emission lines (Si 
IV/O IV], C IV, He II, A1III and UV48 Fe III complex) 
similar to those found bv IR041 ii) the A and B continua 
match very well in all the observed range although the C 
continuum exhibits an excess beyond 2000A; iii) the low 
ionization lines, C III] and Mg II, match reasonably well 
with the same normalization factor as the continua; and 
iv) the EWs of the low ionization lines C III] and Mg II 
are larger in C than in A, and the EW of the C IV line 
is larger in C than in B. 

3.2. Comparison with spectra of other epochs 

Seven spectroscopic observations were made on be- 
tween 2003, May 31 and 2004, March 26 (See Tabled). 
The differences between the WHT spectra correspond- 
ing to componen ts A and B are remarkably similar to 
those reported bv lR04l and lQguri et all 1)200411 from Keck 
observations taken on 2003 May 31. A direct compari- 
son, not shown here (See figures and , among our 
spectra and the Keck 2003 data shows a very good global 
agreement for components A, B, and C. In the case of 
component A, some differences are found in C IV (the 
EW is larger in the WHT spectrum) and in He II (the 
blueward asymmetry is stronger in the Keck spectrum). 




Rest Wavelength (I) 

Fig. 2.— SDSS 1004+4112 A/B spectra taken on 26 March 2004 
at the ARC telescope. 



A direct comparison (not shown here) of the A and B 
ARC spectra based on the normalization of the continua 
shows that during the epoch 2003 November 21 to 2003 
December 22 the blue excess disappears in the line of 
highest ionization (He II) and severely declines in CI V 
but it is still present in the lower ionization Al III and 
Fe III UV48 lines. The WHT observations show that the 
blue enhancements in component A again reached (on 
2004 January 19) a level comparable to that of the Keck 
observations (2003 May 31). In Figure[2]we see that on 
2004 March 26 the system shows similar features. 

3.3. B-A and C-A magnitude difference curves 



INTEGRAL Field Spectroscopy of SDSS 1004+4112 
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Fig. 3. — B — A magnitude difference curve. Triangles: emission 
lines. Circles: co ntinua. Blue squares: photometric data from 
llnada et alJl2003TV Red: WHT-INTEGRAL. Black: Kcck-LRIS. 
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- A magnitude difference curve. Triangles: emission 
continua. Blue squares: photometric data from 
Red: WHT-INTEGRAL. Black: Keck-LRIS. 



We estimate the flux ratio independently for the 
continuum and the emission lines. For the high ioniza- 
tion lines we have integrated only the blue half-side to 
enhance the asymmetries. To obtain the emission line 
flux we subtract a linear interpolation of the continuum 
from the regions adjacent to the lines. We then com- 
puted the continuum flux over intervals of 200 A. The 
resulting magnitude differences, Arris a = tub — rriA and 
AmcA — fnc ~ fnA, for both the continuum and the 
emission lines are presented in Figures 01 and 0J 

In Figure we present the WHT and Keck data. The 
agreement between both observations is quite remark- 
able. The small offset (^0.05 mag) between the Keck 
and WHT continuum difference curves strongly con- 
strains microlensing-induced continuum variability be- 
tween both epochs (2003 May 31 and 2004 January 19). 
There is also good agreement (better than 0.05 mag with 
respect to the WHT continuum curve) with the photo- 
metric data by 103. This is reasonable also in terms of in- 



iH l2002fl to justify chro- 



trinsic variability if we take into account that the yearly 
intrinsic variation of this object could be of about 0.2 
mag (according to the preliminary results of the pho- 
tometric monitoring our group is performing), and that 
the predicted time delay between components A and B 
is about 1 month. 

The B — A differences for the low ionization lines (C 
III], Mg II) match (within errors) the continuum diffe- 
rence curve (there is some marginal evidence of a ne- 
gative offset for the Mg II). This result does not sup- 
port the existence of microlensing in the continuum. The 
blue half-sides of the high ionization lines are above the 
continuum with an offset of about 1 mag for the Si IV, 
C IV, Al III, and the UV48 Fe III complex. The blue 
half-side of the He II line shows a larger offset of about 
3 mag respect to the continuum. Finally, there is no ev- 
idence of wavelength dependence of the continuum diffc- 
rences curve th at have been invoked in several cases (see 
iWucknitz et al 



2003; Mc 

matic microlensing. This also excludes the possibility of 
differential dust extinction. 

The C — A difference magnitude curves are shown in 
Figure 0] The agreement between the WHT difference 
curve for the continuum and the photometric data by 
103 is also remarkably good and does not support the 
existence of noticeable microlensing variability during 
the period between both observations. The Keck spec- 
trum needed an offset of 0.8 mag to be matched to that of 
the 103 data. In Figure0]the high ionization lines appear 
above the continuum but the low ionization lines do not 
match the continuum curve, which is ~ 0.45 mag above. 
An offset between the line and continuum quotients of 
two components of a gravitational len s has been inter- 
prete d as evidence for microlensing (see lMediavilla et"al1 
120051: iMotta et al.l f2002). If we assume that this hypo- 
thesis can be applied to the low ionization BEL and take 
into account that microlensing in the component A is 
for the same reason unlikely, we should accept that the 
continuum of component C is demagnified by microlens- 
ing. On the other hand, the continuum difference curve 
is not flat but shows a relative decrease in C with re- 
spect to A for shorter wavelengths (the color difference 
between components A and C has already been reported 
during the identification of the system as a multiple im- 
aged QSO of large separation, see 103). This can be easily 
interpreted as dust extinction of component C (see e.g. 
Falco et al. 1999 & Muhoz et al. 2004). An interest- 
ing alternative possibility to explain the color differences 
around 2100A and longer (See Figure^) is UV Fell emi- 
ssion microlensed (magnified) in the C image and not in 
the others. The empirical iron emission template given 
by Vestergaard & Wilkes (2001) matches qualitatively 
the observed variation with wavelength. Both possibili- 
ties can be discriminated with future observations, since 
microlensing in the iron lines would show time variability 
but extinction would not. 

4. B-A VARIABILITY FOR THE HEII AND CIII] EMISSION 

LINES 

In Figure \5\ we present the B — A values for the C 
III] and the Hell (blue half-side) emission lines corres- 
ponding to the seven epochs with spectra available (see 
Section 3.2). The B — A value for the low ionization C 
III] line is almost constant (matching the B — A value for 
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the continuum; see previous section), whereas the He II 
emission line exhibits strong variability from a maximum 
during the first epoch, a posterior decay, and a final rise. 
It is very noticeable that the minimum of B — A for the 
He II coincides with the baseline defined by the CIII] . 



-1 




100 200 300 

MJD - 2452790.5 (cays) 

Fig. 5. — Magnitude differences between component A and B for 
CIII] (triangles) and He II lines (squares) 



5. DISCUSSION AND CONCLUSIONS 

The interpretation of the apparent recurrence of the 
blueward enhancement of the high ionization lines is chal- 
lenging. If the enhancements were produced by intrin- 
sic variability, they should appear with some delay in 
all components. However, observations do not reveal 
changes in the spectrum B that could be related to the 
blue wing asymmetries detected in A. On the other hand, 
the lack of continuum variability that may induce the 
change in the emission lines is also against this hypothe- 
sis. 

Alternatively, the rising and fading of the blue wing 
enhancements can be explained by the relative displace- 
ment of a caustic network with respect to the BLR. How- 
ever the displacement of the caustic across a kinemati- 
cally structured BLR should induce velocity changes in 
the high ionization lines l|R04f) rather than recurrences 
like the one observed. Another inconvenience of the mi- 
crolensing hypothesis is the very short time-scale for emi- 
ssion line variability. According to the observations, the 
He II blue wing excess was not present in the ARC 2003 



December 22 spectrum but appears in the WHT spec- 
trum taken on 2004 January 19 (See Figure [SJ implying 
a rise time of less than one month, wh ich re duces by a 
factor of ~6 the upper limit inferred bv IR04I for the size 
of the amplified region under the hypothesis of microlens- 
ing. Richards et al. (2004b) estimated that the impact 
of microlensing in the A continuum was no more than 
20% during the time spanned by their obser vations. Our 
compa rison (see previous section) among the llnada et alJ 
( 2003), Keck, and WHT data indicates no traces of B-A 
continuum variability during the period between 2003 
May 2 and 2004 January 19. On one hand, this im- 
plies that the continuum microlensing time scale (if any) , 
should be much larger than the one associated with high 
ionization line variability. On the other hand, (a) the 
agreement among the B — A data corresponding to the 
low ionization lines and the continuum and (b) the ap- 
parent absence of microlensing chromaticity do not sup- 
port this hypothesis. Thus, it seems that the process 
causing the variability of the high ioniza tion lines does 
not af fect t he co ntinuum. According to iLewis fc Ibatal 
(2004) and R04, microlensing by a caustic network can 
cause uncorrelated variability between the high ioniza- 
tion lines and the continuum. However, it is difficult to 
explain the absence of significant continuum microlens- 
ing (A(B — A) < 0.05 mag) when the high ionization 
lines have undergone a complete cycle that in the case of 
the blue wing of He II implies A(B — A) differences of 
several magnitudes. 

The decoupling between continuum and emission lines 
can be explained by supposing that the emission line re- 
gion is formed by clouds of ionized gas moving at high 
velocity and/or considering a special geometry (biconi- 
cal) for the BLR. This would favor the appearance of 
asymmetrical features and might explain the uncorre- 
lated variability between emission lines and continuum. 
Even the possibility of microlensed transient events as 
SNs or shocks deserves consideration. We will explore 
these possibilities in a future work. 

We are grateful to Gordon T. Richards for providing 
previous observations and useful comments. We also 
wish to acknowledge the referee, Ian Browne for valuable 
comments and ideas. This work was supported by the 
European Community's Sixth Framework Marie Curie 
Research Training Network Programme, Contract No. 
MRTN-CT-2004-505183 "ANGLES", and by the Minis- 
terio de Education y Ciencia of Spain with the grants 
AYA2004-08243-C03-01 and AYA2004-08243-C03-03. 
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